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ABSTRACT 

We investigate the enhancement of star formation efficiency in galaxy interactions and mergers, by numerical simulations of several 
hundred galaxy collisions. All morphological types along the Hubble sequence are considered in the initial conditions of the two 
colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Diff'erent types of orbits are simulated, direct and ret- 
rograde, according to the initial relative energy and impact parameter, and the resulting star formation history is compared to that 
occuring in the two galaxies when they are isolated. Our principal results are: (1) retrograde encounters have a larger star formation 
efficiency (SFE) than direct encounters, (2) the amount of gas available in the galaxy is not the main parameter governing the SEE 
in the burst phase, (3) there is an anticorrelation between the amplitude of the star forming burst and the tidal forces exerted per 
unit of time, which is due to the large amount of gas dragged outside the galaxy by tidal tails in strong interactions; (4) globally, the 
Kennicutt- Schmidt law is retrieved statistically for isolated galaxies, interacting pairs and mergers; (5) the enhanced star formation 
is essentially occurring in nuclear starbursts, triggered by inward gas flows driven by non-axisymmetries in the galaxy disks. Direct 
encounters develop more pronounced asymmetries than retrograde ones. Based on these statistical results, we derive general laws for 
the enhancement of star formation in galaxy interactions and mergers, as a function of the main parameters of the encounter. 



1. Introduction 

Theories about formation and evolution of galaxies have under- 
gone drastic changes in the last decade. Impressive progress in 
the observations of galaxies at high and interme diate redshift 
now put constraints directly on galaxy evolution ( Steidel et al. 
79991 |Le Fevre et aLl[2000l [Cimatti et aL]|2QQ2l ]FTores et al. 



2Q04 |). Th e star format ion history can be estimated (|Madau et 
al. II 1996: Li lly et al.|1 996), and the confrontation with numeri- 
cal simulations help to discover the actual mechanisms for mass 
assembly and transformation of gas into stars ( |Kauffmann &| 
|Charlot|[T998l|Somerville et al.|200T] ). Models are do ne on one 
hand through cosmological N-body simulations (e.g. |Springel 
& Hemquist| |2003[ jSommer-Larsen et al.||2003[ |Governato et 
aL|j2Q04) but in general they lack spatial resolution on galaxy 
scales to resolve star-formation related processes in a realistic 
way; on the other hand, simulations are then completed by semi- 
analytic computations, where galaxy evolution is treated with 
recipes taking into account all important processes, including hi- 
erarchical merging and sta r formation ( jPevriendt & Guiderdoni 
|2000HKauffmann & Haeh nelt 2000|. 

In both cases, it is necessary to study in more details the star 
formation efficiency on galaxy scale, to calibrate the physical 
recipes: star formation rate, and feedback processes have been 
computed on isolated and merging galaxies, and a large vari- 
etyof results have been obtained (Mihos & Hernquist 1994a|b| 



19961 [Sprin gel 2000 ; fTissera et alJ ]2QQ2; Meza et al. 2003; 
Kapferer et al.,2005y|Cox et al.|200 6i All studies demonstrated 
that galaxy collisions trigger star formation, and that the initial 
disk stability is the main parameter to influence the star forma- 
tion sequence: late-type galaxies without bulge are more prone 



to violent bar instability during an encounter, which drives the 
internal gas towar ds the ga laxy center to trigger a nuclear star- 
burst (Mihos & Hernquist 1 1996| ). However, the availability of 
gas in the interacting galaxies is also one of the most determin- 
ing parameter, as well as the adopted gas physics ( |Cox et aL 
[2004). 
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The amount of triggered star formation due to the interac- 
tion is a fundamental parameter required by semi-analytical 
simulations, that strongly depends on models adopted for star 
formation rate and feedback. We propose here to investigate 
this problem statistically, simulating all galaxy types along the 
Hubble sequence, and exploring all physical parameters for the 
interactions, mass ratios, geometrical orbital parameters, etc.. 
in order to get an insight into the global phenomena, and their 
range of variations. In this first paper, the numerical method 
and galaxy models are presented, and star formation evolution 
involving giant-like galaxies, are presented and discussed. A 
wider range of interactions and mergers involving the whole 
mass spectrum of galaxies will be considered in a companion 
paper. Galaxies are initiated at z=0, while simulations involving 
galaxies at any redshift will be studied in a following work. 



The scheme of the paper is as following: in Section [2] the 
adopted galaxy mod els (Section |2.1| ) and the initial orbital 
parameters (Section |2.2|) are described; the numerical code 
is described in Section [3J in Section |4j the main results are 
presented, in particular those related to the evolution of the star 
formation rate during the diflferent phases of interaction, its 
dependence on the total gas amount in the galaxies and on the 
main orbital parameters of the encounter. Finally, in Section |5] 
we analyze and discuss the interactions- starbursts connection 
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Table 1. Galaxy parameters. The bulge and the halo are mod- 
elled as Plummer spheres, with characteristic masses Mb and 
Mh and characteristic radii rg and vh- and Mg represent the 
masses of the stellar and gaseous disks, whose vertical and radial 
scale lengths are given, respectively, by and a*, and hg and Ug. 





gEO 


gSa 


gSb 


gSd 


Mb [2.3 X 10^ Mq] 


70 


10 


5 





Mh [2.3 X lO'^Mo] 


30 


50 


75 


75 


M, [2.3 X 10^ Mo] 





40 


20 


25 


MglM, 





0.1 


0.2 


0.3 


vb [kpc] 


4 


2 


1 




vh [kpc] 


7 


10 


12 


15 


a* [kpc] 




4 


5 


6 


K [kpc] 




0.5 


0.5 


0.5 


Ug [kpc] 




5 


6 


7 


hg [kpc] 




0.2 


0.2 


0.2 



and in Section|6]the main conclusions of this work are drawn. 

2. Initial conditions 

2.1. Galaxy models: moving along the Hubble sequence 

Aiming at exploiting a large set of interactions, involving galax- 
ies of all morphologies from ellipticals to late-type spirals, the 
galaxy models adopted consist in a spherical dark matter halo, 
containing or not a stellar and a gaseous disk and, optionally, a 
central bulge. 

For each galaxy typ^ the halo and the optio nal bulge are mod- 
eled as a Plummer sphere ( [Binney & Trem aine 1987, pag.42), 
with characteristic masses Mb and Mh and characteristic radii 
rB and rn- Their densities are given, respectively, by: 

5/2 

1 + (1) 



pH{r) 
and 
pB{r) 



( 3Mh 
[ Anrn^ 



3Mb 



-5/2 



(2) 



The stellar and gaseous disks follow a IMiyamoto-Nagai density 
profile ( Binney & Tremaine, 1987, pag.44): 



An 



al + (a. + V?T7^7)T' (z2 + h^^f^ 



(3) 



h_lMg 
An 



agR^ -h (a^ -h 3 + hg^) {ag -h + hg^^ 
a] + (ag + ^z^ + /z/) (z^ + hg'^f^ 



(4) 



^ Hereafter we will adopt the following nomenclature for the diff'erent 
morphological types: gEO for giant-like ellipticals, gSa for giant-like Sa 
spirals, gSb for giant-like Sbc spirals and gSd for giant-like Sd spirals. 



Fig. 1. Hubble sequence for our galaxy models. From left to 
right, projection on the x-z plane of a gEO, gSa, gSb and gSd 
galaxy. All the diff'erent components (gas, stars and dark matter) 
are plotted. Dark shading represents densest regions. Each frame 
is 20 kpc X 20 kpc in size. 




Fig. 2. Initial rotation curves for the gSa (solid line), gSb (dashed 
line) and gSd (dot-dashed line) galaxy. 

Table 2. Particle numbers for each galactic composant 





gEO 


gSa 


gSb 


gSd 


N 

^ ^ gas 
star 

Ndm 


80000 
40000 


20000 
60000 
40000 


40000 
40000 
40000 


60000 
20000 
40000 



with masses and Mg and vertical and radial scale lengths 
given, respectively, by K and a^, and kg and ag. Moving 
along the Hubble sequence, from giant-like ellipticals (gEO) 
to giant-like Sd spiral (gSd), the mass of the central spheroid 
varies from Mb = 1.6 x lO^^M© for a gEO to Mb = for a gSd, 
while the gas mass Mg, absent in the case of a gEO, increases 
from 9.2 x IG^Mq in a gSa to 1.7 x IO^^Mq for a gSd (see Table 
[T]for a complete list of all the parameters in Eqs[T]to|4]and Fig[T] 
for a representation of our galaxy sequence). 
The initial rotation curves for spiral galaxies are shown in Fig|2l 



In accordance with observations ( Roberts & Hayne^|1994| ), 



for a fixed distance from the galaxy center, the value of Vrot is 
higher for early type systems than for late type ones, indicating 
a decrease in the enclosed mass content. Aiming at investigating 
interactions between giant like galaxies, the mass ratios of the 
interacting systems is always of order unity. In the next future, 
we will exploit collisions involving galaxies of the whole mass 
spectrum. 

Intending to obtain hundreds of simulations, each galaxy is 
made up of 120000 particles, distributed among gas, stars and 
dark matter, depending on the morphological type (see Table [2]). 

To ini tialize particle ve locities, we adopted the method de- 
scribed in Hernquist ( |1993 ). 
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Table 3. Galaxies orbital parameters 



tree method ( Barnes & Hut|1986" ) and gas evolutio n is followed 



id 


[kpc] 


[kpc] 


[lO^kms-i] 


[lO^km^s-^] 


spin'^ 


Oldir 


100. 


8.0 


7.07 


0.0 


up 


Olret 


100. 


8.0 


7.07 


0.0 


down 


02dir 


100. 


8.0 


7.42 


2.5 


up 


02ret 


100. 


8.0 


7.42 


2.5 


down 


03dir 


100. 


8.0 


7.74 


5.0 


up 


03ret 


100. 


8.0 


7.74 


5.0 


down 


04dir 


100. 


8.0 


8.94 


15.0 


up 


04ret 


100. 


8.0 


8.94 


15.0 


down 


05dir 


100. 


16.0 


5.00 


0.0 


up 


05ret 


100. 


16.0 


5.00 


0.0 


down 


06dir 


100. 


16.0 


5.48 


2.5 


up 


06ret 


100. 


16.0 


5.48 


2.5 


down 


OVdir 


100. 


16.0 


5.92 


5.0 


up 


07ret 


100. 


16.0 


5.92 


5.0 


down 


08dir 


100. 


16.0 


7.42 


15.0 


up 


08ret 


100. 


16.0 


7.42 


15.0 


down 


09dir 


100. 


24.0 


4.08 


0.0 


up 


09ret 


100. 


24.0 


4.08 


0.0 


down 


lOdir 


100. 


24.0 


4.65 


2.5 


up 


lOret 


100. 


24.0 


4.65 


2.5 


down 


lldir 


100. 


24.0 


5.16 


5.0 


up 


llret 


100. 


24.0 


5.16 


5.0 


down 


12dir 


100. 


24.0 


6.83 


15.0 


up 


12ret 


100. 


24.0 


6.83 


15.0 


down 



^ For two equal point masses of mass m = 2.3 X 10 M©. 
^ It is the total energy of the relative motion, i.e. 

E - v^/2 - G(mi + m2)/r. 
^ Orbital spin, if parallel (up) or antiparallel (down) to the galaxies 

spin. 



2.2. Orbital parameters 

Aiming at exploiting a vast range of orbital parameters, for each 
couple of interacting galaxies we performed 24 different simu- 
lations, varying the galaxies orbital initial conditions, in order to 
have (for the ideal Keplerian orbit of two equal point masses of 
mass m = 2.3 X 10^ ^M©) the first pericenter separation r^^^ =8, 
16, and 24 kpc. For each of these separations, we varied the rel- 
ative velocities at pericenter, in order to have one parabolic and 
three hyperbolic orbits of different energy. Finally, for each of 
the selected orbits, we changed the sign of the orbital angular 
momentum in order to study both direct and retrograde encoun- 
ters. 

Combining each orbital configuration with all possible mor- 
phologies for the interacting pair of galaxies, we obtained a total 
sample of 240 interactions. 

In Table [3] the initial distance r/„/ and the pericenter distance 
rp between the galaxies center-of-mass are listed, together with 
their relative velocity Vp at pericenter and the orbital energy E, 
for all the simulated encounter^ 

The results described in this paper refer only to planar encoun- 
ters (galactic disks lying in the orbital plane). 



3. Numerical method 

To model galaxy evolution, we employed a Tree-SPH code, in 
which gravitational forces are calculated using a hierarchical 



^ The values refer to the ideal Keplerian orbit of two equal point 
masses of mass m - 2.3 x lO^^M©. 



by means of smoothed particle hydrodynamics ( |Lucy| \T9T7\ 
Gingold & IMonaghan 1982). Gravitational forces are calculated 
using a tolerance parameter 6 = 0.7 and including terms up to the 
quadrupole order in the multiple expansion. A Plummer poten- 
tial is used to soften gravitational forces, with constant softening 
lengths for different species of particles. In the simulations de- 
scribed here, if not explicitely otherwise indicated, we assume 
6 = 280pc. 

SPH (Smoothed particle hydrodynamics) is a Lagrangian tech- 
nique in which the gas is partitioned in fluid elements repre- 
sented by particles, which obey equations of motion similar to 
the collisionless component, but with additional terms describ- 
ing pressure gradients, viscous forces and radiative effects in gas. 
To capture shocks, a conventional form of the artificial viscosity 



is used, with parameters a = 0.5 and JS = l.O (Hernquist & Katz 



1989). To describe different spatial dynamical range, SPH par- 
ticles have individual smoothing lengths hi, calculated in such 
a way that a constant number of neighbors is contained within 
2hi. All the simulations have been performed using a number of 
neighbors A^^ ~ 15. The gas is modeled as isothermal, with a 
temperature Tgas = 10"^ K. Because of the short cooling time of 
disk gas, fluctuations in the gas temperature are quickly radiated 
away, so that simulations employing an isothermal equation of 



state differ little from more realistic ones ( jlMihos & Hernquist 
1996l|Naabetal.|20 06D. 

The equations of motion are integrated using a leapfrog algo- 
rithm with a fixed time step A^ = 5 x lO^yr. 

3.1. Star Formation and continuous stellar mass loss 

Including star formation in gas dynamics is not a trivial task and 
a lot of different recipes and numerical methods can be adopted 
([Katz 1992; Steinmetz & JMiiller 1994; Sp ringel 2000,; Springel] 
& Her nquist 20T)3f|Cox et al.||2Q Q6), in order to model the star 
formation rate on one hand, and taking into account the effects 
that this star formation has on the surrounding, on the other hand. 



As in |]V[ihos & Hernquist ( |1994b| ), we parametrized the star 
formation eflficiency for a SPH particle as 



= CXi 



1/2 



(5) 



gas 



with the constant C chosen such that the isolated disk galaxies 
form stars at an average rate of between 1 and l.SM^yr'^. 
The choice of the parametrization in Eq|5] is consistent with 
the observational evidence that on global scales the SFR in 
disk galaxies is well represented by a Schmidt law of the form 
being Ugas and Us fr disk-averaged surface den- 



■'SFR 



'gas 



sities, with the best fitting slope about 1.4 (s ee Kennicutt 
T998a| but also|Wong & Blitz|2002t|Boissier et al. 2003 ; Ga£& 
Solomon 2004). Interestingly, this relation seems to apply, with 
a simil ar slope, also to local scales, as shown in , Kennicutt et al.| 
( [2005] ) for M51. 

Once the SFR recipe is defined, we apply it to SPH particles, us- 
ing the hybrid method described in IMihos & Hernquist (1994b]): 
it consists in representing each gas particle with two mass val- 
ues, one refering to its gravitational mass M/, whose value stays 
unchanged during the whole simulation, and the other describ- 
ing the gas content of the particle Mt^gas, whose value changes in 
time, according to Eq|5] Gravitational forces are always evalu- 
ated on the gravitational mass M/, while hydrodynamical quanti- 
ties, in turn, uses the time-varying Mi^gas- Only if the gas fraction 
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present into the hybrid particles drops below the 5% of the initial 
gas content, the hybrid particle is totally converted into a star- 
like particle and the little amount of gas material still present is 
spread into the neighbors. 

We also followed the method described in IMihos &1 
Hernquist] ( | 1 994 b ) for including the effects of star formation into 
the ISM. The method is fully described in the above cited paper, 
and here we briefly recall only the main features. A Miller- Scalo 
stellar mass function is adopted, and we evaluated the fraction 
of stars with masses > SM©, assuming that they instantaneously 
become supernovae, living behind remnants of 1.4M0 and re- 
leasing their mass on the surrounding ISM. The mass released 
enriches also the metallicity of the surrounding gas. This is done 
assuming a yield y = Mret/M^=0.02, Mret being the total mass 
of all reprocessed metals and the total mass in stars. For each 
gas particle, mass and metals return is applied to the / - th neigh- 
bor gas particle, using a weight w/ based on the smoothing ker- 
nel. 

The energy injection in the ISM from SNe explosions is treated 
assuming that only a fraction 6kin of Esn = 10^^ erg goes into ki- 
netic energy, by applying a radial kick to velocities of neighbor 
gas particles; thus, for each SNe explosion, the i-th neighboring 
gas particle receives a velocity impulse directed radially away 
from the "donor", with a magnitude 



Avi = 



I 2 Wi Ekin 



SN 



Mi 



1/2 



(6) 



Wi being, once again, the weighting based on the smoothing 
kernel and Mi the mass of the receiver. 

Evidently, this method has a certain number of free parame- 
ters, that clearly influence not only the star formation evolution, 
but, more generally, the global galaxy dynamics. Before moving 
on to perform simulations of galaxy encounters, we run a 
set of simulations of isolated galaxies, in order to check the 
dependence of the results on the gravitational smoothing length 
6, and on the fraction ekin of kinetic energy that is released, 
via SNe explosions, to the surrounding gas particles. For 
example, for the choice of the gravitational smoothing length, 
we performed 21 simulations of galaxies gSa, gSb and gSd, with 
6 = 70,140,210,280,350,420 and 490 pc, respectively. The 
results showed that a too small value of e causes a conspicuous 
heating of the stellar disk, for all the morphological galactic 
types. An 6 = 280/7C, which corresponds (for gSa and gSb 
galaxies) to the average distance among stars in the disk seemed 
to represent a compromise between accuracy in the gravitational 
forces evaluation and not excessive relaxation effects. Also for 
the choice of the Ekin parameter, a set of numerical simulations 
of isolated galaxies was perfo rmed, and the results found were 
in good agreement to those of |Mihos & Hernquist| ( |1994b| ), i.e. 
it was checked that if the total amount of kinetic energy received 
by a gas particle, due to the contribution from all its neighbors, 
is < lkms"\ this prevents the gaseous disk from a rapid growth 
of the vertical thickness, giving good results in terms of gas and 
star formation morphology. 

If star formation acts in consuming the gas of a galaxy, an 
exhaustive modelling of galactic evolution has to take into ac- 
count also the competing process of stellar mass-loss. Indeed, as 
both observational studies and stellar evolutionary models pro- 
pose, it can play an important role, since the gas mass fraction 
restituted by stars may reach some 45% over the Hubble time, 
when integrated over the stellar mass spectrum. 



We included continuous stel lar mass loss in the mo del adopting 
the formula given in Eq.2 of Jungwi ert et aLl pOOl), and apply- 
ing it only to stellar populations formed into hybrid particles, 
i.e., at each time step, an amount 



Mi,s{t) 



[Mi - Mi^gasit)) At Co 
t - hirth + ^0 



(7) 



of the stellar mass in the population is lost by evolutionary ef- 
fects, going to enrich the gas content Mi^gas of the hybrid particle. 
In the formula above, tMrth represents the birth ti me of the popu- 
lation, To = ASlMyr and cq = 5.47 x 10"^ (see [Jungwiert et al. 
2001, for details). 

4. Results 

In this Section, the main results of our study will be presented 
and discussed. This first paper being devoted to the study of 
star formation in interacting galaxy pairs, we will exclude 
from our analysis 24 simulations from the total sample of 240, 
these simulations involving only gEO-gEO encounters ("dry" 
interactions). 

After presenting an image gallery of galaxy interactions and 
discussing the ma in m orphological features arising during 
encounters (Section [4~T] ), we will move to examine the evolution 
of the star formation. To distinguish which contribution to the 
global SFR comes from the undisturbed galaxies and which 
from the interaction, we will first analyze the evolution of the 
star formation rate for the isolated galaxies (Section |4.2| ), then 
we will proceed to describe the most salient features of the SFR 
evolution for the interacting pairs (Section |4.3| ). The evolution 
of the star formation efficien cy (S FE) too will be discussed 
(Section [4.4| ), while in Section [43] we will describe where star 
formation regions are located. 

The next step will consist in trying to understand which 
parameters determine the different SFR obtained soon after 
the pericenter passage and in the merging phase, for different 
galaxy pairs. This will lead us to study in detail the dependence 
of the SFR on: total gas content Mgas in the galaxy just before 
the burst phase (Section 4.6), galaxies separation at pericenter 
passage Rp, galaxies relative velocity at pericenter passage V^, 
characteristic time of the encounter tenc^ tidal effects of the 
encounter, quantified by means of a suitable 'tidal parameter' 
(Section |4.7| ). The discussion on the dependence o f th e SFR 
on this tida l par ameter will be deepened in Section 4^ while 
in Section 4.9 we will propose a formulation for the SFE 
in the burst phase. Section |4.10| presents a discussion of the 
evolution of interacting, merging and post -merg er galaxies into 
the (X^gas.^sFR) plane. Finally in Section 4.11 gas inflow into 
the central galactic regions will be exploiteo 

In the following, we will often refer to galaxy center in our 
analysis. If not otherwise written, we refer to a density- weighted 
center evaluated in the following way: 
for each configuration, and for each galaxy. 



- first we evaluate the density center ( [Casertano & Hut|1985| ) 
of the dark matter particles Cdm initially (t=0) belonging to 
the galaxy; 

- then we evaluate the density center for all the particles 
(gas-Fstars-Fdark matter) initially (t=0) belonging to the 
galaxy, that are at a distance r < 10 kpc from Cdm\ the 
density center so found is used to define the center C of the 
whole galaxy. 
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Fig. 3. Evolution of gas (upper panels), stars (intermediate panels) and dark matter (lower panels) during a direct merger between 
two gSb galaxies (id=09dir in Table [3|. Time is labelled in the upper part of the Figure. Each frame is 50 kpc x 50 kpc in size. 



For evaluating the galaxy velocity, we evaluated the mass- 
weighted velocity of all the particles ( stars + gas + dm) at a dis- 
tance r < lOkpc from C. 

We will also refer sometimes to specific encounters. In this case 
the nomenclature adopted is the following: morphological type 
of the two galaxies in the interaction (gEO, gSa, gSb or gSd), -h 
the encounter identification string (see first column in Tabl^, -l- 
00 (for indicating a planar encounter). For example the nomen- 
clature gE0gSa04ret00 corresponds to a planar interaction be- 
tween an elliptical and a Sa spiral, whose initial orbital parame- 
ters are those corresponding to id=04ret in Tabl^ 

4.1. A gallery of galaxy interactions 

Given the large number of simulations performed, it is not 
possible to describe each case individually. So we will proceed 
to describe some 'fiducial' cases, that retain the mean features 
of all the simulations performed. In FigsjSj [4] and [5] some 
sequences of galaxy mergings and flybys are shown. 

Fig. [3] describes the encounter and successive merger of two 
gSb galaxies, on direct orbit. At the beginning of the simulation 
the two galaxies are separated by a distance of 100 kpc. As 
they start to approach each other, they begin to develop tails, 
populated of both stars and gas particles. The intense tidal field 
exerted during the pericenter passage (t=400 Myr) leads also to 
a transfer of mass between the two systems. In many cases, as 
it will be shown in the next sections, the direct encounter leads 
to a more rapid and dramatic merger, with an expansion of the 
outer parts of the system, which is particularly visible for disk 
galaxies. 

The retrograde merger is shown in Fig.|4] In this case the two 
galaxies are less aff'ected by the tidal interaction. The disk galax- 



isogSa 



isogSb 



isogSd 




1000 2000 3000 
Time (Myr) 



1000 2000 3000 
Time (Myr) 



1000 2000 3000 
Time (Myr) 



Fig. 6. From left to right: Star formation rate of the isolated gSa, 
gSb and gSd galaxies. 



ies develop two transient great spiral arms, after the pericenter 
passage, but retain much longer their initial structure. Contrary 
to the direct case, here no transfer of mass from one galaxy to 
the other takes place and the formation of tidal tails is less obvi- 
ous. As it will be discussed later, this allows most of the initial 
gas mass to stay well confined in the disk of the spiral galaxy, 
furnishing a great reservoir for the intense starburst that takes 
place in the final stage of the merging. Note indeed the high gas 
concentration in the inner central galactic regions in the last two 
snapshots of galaxies in Fig|4] 

Finally, Fig. [5] shows a direct flyby between two gSa spirals. 
Tidal tails and a bridge connecting the two galaxy centers de- 
velop after the pericenter passage at t=170Myr. In this, as in the 
previous Figures, the high density regions in the gas component 
correspond also to the sites of m ost intense star formation, as it 
will be discussed in Section |431 
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Fig. 4. Evolution of gas (upper panels), stars (intermediate panels) and dark matter (lower panels) during a retrograde merger 
between two gSb galaxies (id=09ret in Table jSj. Time is labelled in the upper part of the Figure. Each frame is 50 kpc x 50 kpc in 
size. 
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Fig. 7. Evolution of the gas component (upper panels) and of the 
stellar one (lower panels) for the isolated gSb galaxy. Time is 
labelled in the upper part of the Figure. Each frame is 20 kpc x 
20 kpc in size. 



4.2. Evolution of isolated galaxies 

Since the goal of this paper is to investigate the role interactions 
and mergers play in the star formation process, it is essential to 
study the evolution of the SFR in the isolated galaxies, in order to 
distinguish secular evolution from tidal effects in the interacting 
galaxy sample. This is shown in Fig. [6] In all the three cases, an 
initial burst develops, driven by the compression of the gas into 
density waves. In Fig|7] which shows the time evolution of the 
gaseous and stellar component for the gSb galaxy, spiral arms 
and a bar are clearly formed at t=500Myr. 

After the initial peak, the SFR follows a typical exponen- 
tial profile. This is a general result of simulations of galaxies 
evolved in isolation (see, for example, the discussions in Combes 
2004[ |Harfst et"aLl|2006| ), unless ad hoc recipes for modeling 
feedback from SNe explosions are taken into account ([Springel 



2000| ). To reproduce the average constant star formation rates 
found in observati ons of sp iral galaxies in t he ni iddle of the 
Hubble sequence ((Kennicuttj[1983 ; Kennicutt et al.||1994| ), it is 
indeed necessary to take into account also external gas accre- 
tion, as shown in simulations of g alaxy evolution in a cosmolog- 
ical frame ( |Tissera|2000[ Nagamine et al.|2004j ). 
Averaged on the first 3 Gyr of evolution, the <SFR> varies from 
IMo/yr to 2.5Mo/yr. 

4.3. Star formation in interacting pairs 

The variety of morphological and orbital parameters adopted in 
the simulations reflects into the different star formation histories 
of galaxies in the sample. Some SFR evolution with time are 
shown in Figs[8]and[9j which refer, respectively, to mergings and 
flybys. Simply looking at these few examples, it appears clear 
that interactions can increase star formation in the galaxy pairs 
from low levels (1.5-2 times the isolated case) to values typical 
of starburts galaxies (20-60 times the isolated case). 
In these figures, no distinction is made between the two inter- 
acting galaxies in the sample, i.e. the relative SFR at time t is 
evaluated as the ratio between the SFR of the pair and the to- 
tal SFR of the isolated galaxies, both evaluated at time t. From 
Fig [8] one can extrapolate some of the main features of the star 
formation evolution for interactions leading to mergers: 

- the peak of the SFR usually occurs in the last phases of the 
encounter, when the two galaxies are interpenetrating, while, 
at the first pericenter passage, the increase in the SFR is only 
modest (1 to 4-5 times that of the isolated case). 

- for the same orbit, direct and retrograde encounters can lead 
to different SFR evolution. Obviously, this is not a general 
rule, but often the SFR is greater for retrograde mergers than 
direct ones and in some cases the differences between direct 
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Fig. 5. Evolution of gas (upper panels), stars (intermediate panels) and dark matter (lower panels) during a direct flyby between two 
gSa galaxies (id=04dir in Table [3|. Time is labelled in the upper part of the Figure. Each frame is 50 kpc x 50 kpc in size. 



and retrograde encounters is conspicuous. For example in the 
gEOgSa05retOO encounter, the maximum SFR is 20 times 
greater than of the direct gEOgSa05dirOO encounter (upper 
panels in the Figure). 



For flybys, the peak in SFR occurs usually just after the 
pericenter passage; also in this case the increase in the star 
formation rate can be modest (as for the gSbgSdOTdirOO 
and the gEOgSblOretOO cases in Fig|9]) or substantial (cf. the 
gSagSa04dir00 and gSbgSbOSretOO in the same Figure). On 
average, the SFR increases by a factor 3.4 with respect to the 
isolated galaxy sample. 



Note that, both in Fig[8 
mation rate corresponds to t 



and Fig|9] the peak in the star for- 
le occurence of a nuclear starburst. 



4.3.1 . Integrated star formation 

Here we want to examine the evolution of the integrated star for- 
mation rate (iSFR = ^^^^^^"^ SFR{t)dt^, i.e. the total quantity of 

gas transformed into stars, during the whole duration of the sim- 
ulations. We have seen that different simulations show different 
SFRs, therefore they are more or less efficient, during pericenter 
passages or mergers, in transforming gas into stellar matter. It is 
interesting to check if the ISFR shows the same trend (for ex- 
ample if retrograde mergers are more efficient than direct ones), 
because a priori this is not assured: ISFR can be dominated by an 
extreme starburst-like event for a short duration (several million 
years) or a long-term low enhancement of SFR due to interac- 
tion. 

In Fig[TT] the ISFR, normalized to the isolated case, is shown, 
for all the simulations performed, as a function of the orbit. The 
main features are: 



To show the wide variety of SFR obtained, in Fig.fTOlwe re- 
port the maximum SFR obtained during the 3 Gyr-evolution of 
the pairs, as a function of the different orbits simulated. As be- 
fore, the SFR is relative to the isolated case. In this plot, results 
are shown for the whole sample (i.e. both mergers and flybys). 
It clearly puts in evidence that, for the same initial orbital condi- 
tions, different morphologies of galaxies in the pair lead to dif- 
ferent SFR histories, but also that, fixing the morphology of the 
pair, direct and retrograde encounters usually do not lead to the 
same star formation evolution. For mergers in general retrograde 
encounters are more efficient that direct ones in driving star for- 
mation: indeed, on average, in retrograde mergers the SFR peaks 
to 8.9 times that of the corresponding isolated galaxies, while, 
for direct mergers, the maximum star formation rate (always rel- 
ative to the isolated case) is only 5. 



- Integrated on 3 Gyr of evolution, the total gas mass converted 
into stars can be increased by up to a factor ~ 2, due to galaxy 
interactions. 

- Flybys can be as eflftcient as mergers in transforming gas into 
stars: for example, two gSb galaxies, involved in a retro- 
grade flyby, can have an ISFR comparable to that obtained 
in a merger (cfr. the L-L flyby with orbit id=03 with the L-L 
mergers 05 and 09 in Fig. 
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- In general, retrograde encounters are slighlty more efliicient 
in transforming gas into stars: on average, for retrograde in- 
teractions < ISFR >= 1.24, while < ISFR >= 1.15 for 
direct encounters. 

- The highest ISFR values are obtained for interactions involv- 
ing late type galaxies. 
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Fig. 8. Star formation rate, versus time, for some galaxy merg- 
ers. The SFR is normalized to that of the corresponding isolated 
galaxies. The blue arrows indicate the first pericenter passage 
between the two galaxies, and the red arrows the merger epoch. 
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Fig. 9. Star formation rate, versus time, for some galaxy fly- 
by s. The SFR is normalized to that of the corresponding isolated 
galaxies. The blue arrows indicate the epoch of pericenter pas- 
sage between the two galaxies. 



4.4. Star formation efficiency in interacting pairs 

To describe how much stars are produced by unit gas mass, it 
is useful to study the star formation efficiency (SFE). This pa- 
rameter can be defined as the ratio between the amount of gas 
transformed into stars, at a generic time t, and the available gas 
content at the same time: 



SFE(t) = 



Moas(t) 



(8) 



Several authors indicate that strong ly interacting gal axies are 
also more efl&cient in forming stars ([Solomon & S age 1988; 
Combes et"aLl 119941 [Georgakakis et al.||2000| ), while others 
( Casasola et al. |2004l i found that even if interacting galaxies ap- 
pear more luminous in the infrared, this higher star formation 
rate does not correspond to a diflTerent efl&ciency in star forma- 
tion per unit of calculated H2 mass. 

In Fig[T2]the SFE, relative to that of the corresponding iso- 
lated galaxies, is shown for some of the simulation performed. 
In all the cases, we found that an increase in the SFR is corre- 
lated with an increase in the SFE, as expected from Eq[8] It is 
interesting to note, however, that, for the most intense bursts: 

- usually the SFE peak is higher than the corresponding SFR 
and is also slightly delayed in time 



- the characteristic SFE times are greater than the correspond- 
ing SFR ones, i.e., after the encounter/merger many galaxies 
can show a SFR which is returned at preinteraction levels, 
while the SFE still shows higher values in comparison to the 
isolated counterparts. 

This is evidently a consequence of the adopted SFE defi- 
nition, because of the fact that, during an intense burst of star 
formation, gas in interacting galaxies is more rapidly depleted, 
with respect to the isolated case, thus the total amount of gas 
content in the galaxies is then lower than that of the isolated 
counterpart. 



In Fig. 13 the maximum star formation efficiency obtained 
for each encounter is shown, as a function of the orbit. 



4.5. Where do stars form? 

From the results presented in the previous sections, it is clear 
that interactions and mergers can lead to a conspicuous en- 
hancement both in the SFR and in SFE (up to a factor of 100 
for the SFE relative to the isolated case), but this is not always 
so. The picture that begins to emerge is that interactions, and 
mergers in particular, can lead to intense starburst but this is 
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Fig. 10. Maximum star formation rate, relative to the isolated 
case, as a function of the orbital type, for flybys (empty symbols) 
and mergers (solid symbols). Encounters are grouped into three 
classes, depending on the morphology of the interacting galax- 
ies: interactions between two early type galaxies are shown in 
the upper panel, the medium panel refers to early-late type en- 
counters, the lowest to encounters involving two late-type sys- 
tems. In each panel, error bars represent the standard error of the 
mean. The large error bars found in some cases reflect a large 
dispersion in the data. 
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Fig. 11. Integrated star formation rate, ISFR, relative to the iso- 
lated case, as a function of the orbital type, for flybys (empty 
symbols) and mergers (solid symbols). Encounters are grouped 
into three classes, depending on the morphology of the interact- 
ing galaxies: interactions between two early type galaxies are 
shown in the upper panel, the medium panel refers to early-late 
type encounters, the lowest to encounters involving two late-type 
systems. In each panel, error bars represent the standard error of 
the mean. The large error bars found in some cases reflect a large 
dispersion in the data. 



not always the case, i.e. galaxy interactions are not a sufficient 
condition to convert high gas mass quantities into new stars 
(see [Bergvall et al. 2003 for an observational study of a sample 
of interacting galaxies that result in poor starburst triggers). 
While, up to now, numerical simulations seem to support 
mainly the idea of a frequent occurence of tidally triggered 
central starbursts (| Mihos &Tlernquist 1994a|b| |1996t |Springel| 
'2000| ), observational results draw a more complex scenario. 
Ultra Luminous Infrared Galaxies (ULI RGs) indicate centra lly 
concentrated star formation (but see also [Combes et al.|2006| for 
an example of a moderate redshift ULIRG with a more extended 
CO emission), while a certain number of interacting galaxies 
presents a more extended star formation. In the Antennae 
galaxies ( [Wang et aL]|2004| ), for example, the most intense star 
formation regions are located between the two galaxies, and 
the late-type galaxy NGC 275 in the Arp 140 system, despite 
being close to a merger with its companion NGC 274, does not 
show evidence of an enhanced or centrally concentrated star 
formation, the brightest emissio n from star form ation tracers 
coming from an off-center region ( [CuUen et al.|2 006 ). 



Figs 14 [15 16 17 18 and 19 show a comparison of most 



intense star formation regions and gas maps for three mergers of 
the sample. Far to represent the whole sample, they clearly sug- 
gest that the tidally triggered star formation process can manifest 
in a variety of ways. The retrograde encounter/merger between 
two late type Sb galaxies leads to a strong enhancement in the 



totaj^ SFR (about 20 times that of the isolated case), the site 
of most intense star formation being strongly centrally concen- 
trated. The same encounter, this time involving galaxies in direct 
orbits, leads to the star formation maps shown in Fig 16 in this 
case the central galaxy regions are sites of intense star formation 
in the final phases of the merging process, but, this time, the tidal 
tails and the bridge connecting the two galaxies contribute also 
to the production of new stars, as well as circumnuclear ring^ 
The situation drastically changes in Fig. 18 where the star for- 
mation maps resulting from an encounter between an early-type 
elliptical and a late-type Sd galaxy are shown. In this case, the 
first pericenter passage between the two galaxies is highly dis- 
ruptive for the gSd, which develops two giant tails, that acquire 
a high fraction of gas mass, and that are also site of intense, lo- 
cal, and clumpy star formation. Note that the star forming clump 
in the upper left part of the panels in Fig corresponds to the 
center of the elliptical galaxy, that has acquired part of the gas 
mass of the companion during the interaction. 
On average, the gas mass being spread at great distances from 
the galaxy center, the gas density is obviously lower than the one 
obtained in the above cited encounter. This lead to a total SFR 
which is sensibly lower than the one obtained in the other two 
mergers: only 3 times higher than the isolated case. 



^ i.e. evaluated on the total gas mass of the galaxy. 

^ Note that the presence of an ILR stops the gas inflow. This, together 
with the tidal gas removal, explains why SFRs tend to be lower in direct 
orbits, as we will see later on. 
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Fig. 12. Star formation efficiency (cyan curve), relative to 
the isolated case, as a function of time for three mergers 
(gEOgSaOSretOO, gEOgSbOlretOO, gSbgSb09ret00) and a flyby 
(gSdgSdOSdirOO). The blue arrows indicate the first pericenter 
passage between the two galaxies, and the red arrows the merger 
epoch. For comparison also the relative SFR is shown (grey 
curve) 



4.6. How does the maximum SFR depend on the total 
amount of gas available in the galaxies? 

Before moving on to describe in details the gas and stars 
dynamical evolution during the encounters, in order to try to 
understand the physical processes that lead to an enhancement 
in the SFR and in the ISFR, with respect to the isolated cases, 
we want to understand whether, and to what extent, the SFR 
peaks depend on the amount of gas available in the galaxy. In 
other words, are the diff'erent SFR peak values shown in Figs. [8] 
|9] [TOj [TT] and [12] due to a diff'erent gas content in the galaxy at 
the moment of the burst? 

To investigate this, we evaluated the gas contenj^ (relative to 
that at the beginning of the simulation) just before the burst, i.e. 
for flyby s we evaluated the amount of gas in the galaxy 50 Myr 
before the pericenter passage, while for mergers, the quantity 
is evaluated 50 Myr before the coalescence between the two 
galaxies. In Fig|20 the maximum SFR is plotted, as a function 
of the gas content, expressed in units of the initial galactic gas 




Fig. 13. Maximum star formation efficiency, relative to the iso- 
lated case, as a function of the orbital type, for flybys (empty 
symbols) and mergers (solid symbols). Encounters are grouped 
into three classes, depending on the morphology of the interact- 
ing galaxies: interactions between two early type galaxies are 
shown in the upper panel, the medium panel refers to early-late 
type encounters, the lowest to encounters involving two late-type 
systems. In each panel, error bars represent the standard error of 
the mean. The large error bars found in some cases reflect a large 
dispersion in the data. 



mass, for all the encounters. Flybys and mergers locate in two 
difl'erent regions of the plots: for flybys, the pericenter passage 
occurs in the initial phases of evolution, when the gas content of 
the disk galaxies is similar to the initial ones, while the mergers 
are located in a region of space where the amount of gas mass 
is a fraction between 0.2 and 0.4 of the initial gas mass for 
elliptical- spiral encounters, while it is in a somewhat larger 
range (0.1-0.6) for spiral-spiral ones. 

In both cases (flybys and mergers) it is evident that one can have 
large difl'erences in the maximum SFR, even when the amount 
of fuel available is nearly the same: this is the case of flybys 
involving a gEO and a gSb in direct orbits, for example, where 
one can have a factor of 4 in the maximum SFRs, even if the gas 
mass is comparable. Even more striking is the case of mergers 
where, with an available gas mass which is about 33% of the 
initial value, there are encounters which lead to a maximum 
SFR 60 times the isolated case, while others enhance the SFR 
only of a factor of 5. 

In conclusion. Fig |2Q| clearly shows that the quantity of fuel 
still available in the galaxy is not the main parameter which 
influences the SFR in the burst phase. 



4.7. How does the maximum SFR depend on Rp,Vp, tenc ^nd 
on the effect of tidal forces? 



^ The gas content considered is the total amount of gas, i.e. the gas 
can be present both in the tails or in the main body of the galaxy. 



In the previous Section, we have found that the maximum SFR 
in interacting and merging galaxies does not depend mainly on 
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Fig. 14. Star forming regions during the final stage of an 
encounter involving two gSb galaxies in retrograde orbits 
(id=09ret in Table[3|. Only hybrid particles with a star formation 
efficiency (averaged on 50 Myr) greater than 0.005 are shown. 
Each frame is 20 kpc in length. Snapshots are shown every 50 
Myr. Blue-green colors correspond to the highest densities re- 
gions. See Fig[T5]for a comparison with gas maps. 



Fig. 16. Star forming regions during the final stage of an en- 
counter involving two gSb galaxies in direct orbits (id=09dir in 
Table [3|. Only hybrid particles with a star formation efficiency 
(averaged on 50 Myr) greater than 0.005 are shown. Each frame 
is 20 kpc in length. Snapshots are shown every 50 Myr. Blue- 
green colors correspond to the highest densities regions. See 
Fig [17] for a comparison with gas maps. 





- ^ 




, 1 1 1 1 1 , , , 1 


: ' 1 ' " 1 ' — ^ 


H ^ ^ 


^ ^ 


1 


: ' 1 ' " 1 ' — ^ 


H ^ ^ 


^ — ^ — ^ 
* 


• 

' 1 ' 



V 


} 






: ' 1 ' " 1 ' — ^ 

^ ^ - 


H ^ ^ 


^ — ^ — ^ 


^ — ^ — ^ 


: ' 1 ' " 1 ' — ^ 


H ^ ^ 

- - 


^ — ^ — ^ 


^ — ^ — ^ 



Fig. 15. Gas maps during the final stage of an encounter involv- 
ing two gSb galaxies in retrograde orbits (id=09ret in Table [3|. 
Each frame is 20 kpc in length. Snapshots are shown every 50 
Myr. See Fig[T4]for a comparison with star forming regions. 



Fig. 17. Gas maps during the final stage of an encounter involv- 
ing two gSb galaxies in direct orbits (id=09dir in Table |3]). Each 
frame is 20 kpc in length. Snapshots are shown every 50 Myr. 
See Fig 16 for a comparison with star forming regions. 



the total amount of gas mass available in the main body of the 
system and in the tails. A first indication of the fact that the 
encounter geometry aff'ects the enhancement of star formation 
has been found in Section |4.3[ comparing retrograde and di- 



rect encounters: we have seen that usually retrograde interac- 
tions/mergers are more efficient than direct ones in converting 
gas to stars. But it is still clear that other parameters must play 
a role in the evolution of the SFR in interacting galaxies: even 
considering only retrograde encounters between galaxies with 
the same morphology, in some cases the increase in the SFR rel- 
ative to the isolated evolution is only modest, while other events 
lead to a SFR 60 times greater than that of the isolated case. 
So, the aim of this Section is to explore the dependence of the 
maximum SFR (occurring in the coalescence phase for mergers 



and soon after the pericenter passage for flybys) on other param- 
eters, such as: 



- the distance at the first pericenter passage Rp ; 



the relative velocity Vp at 
the characteristic time of the encounter, L, 



Rp and Vp have been evaluated using the definitions of galaxy cen- 



ters and velocities given at the beginning of Sect|4] 
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Fig. 18. Star forming regions during the final stage of an en- 
counter involving a gEO and a gSd galaxy in direct orbits 
(id=03dir in Table[3|. Only hybrid particles with a star formation 
efficiency (averaged on 50 Myr) greater than 0.005 are shown. 
Each frame is 100 kpc in length. Snapshots are shown every 50 
Myr. Blue-green colors correspond to the highest densities re- 
gions. See Fig[T9]for a comparison with gas maps. 
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Fig. 19. Gas maps during the final stage of an encounter involv- 
ing a gEO and a gSd galaxy in direct orbits (id=03dir in Table[3]). 
Each frame is 100 kpc in length. Snapshots are shown every 50 
Myr. See Fig 18 for a comparison with star forming regions. 
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Fig. 20. Maximum star formation rate (relative to the isolated 
case), as a function of the gas content (rescaled to its initial 
value), for flybys (empty symbols) and mergers (solid symbols). 
Diff'erent symbols correspond to diff'erent morphologies of inter- 
acting galaxies, as explained in the legend. Error bars represent 
the standard error of the mean. 



to the interactio ns with a compan ion galaxy whose mass is 
^comp (see also Boumaud et al.|2005^ ). As a measure of the 
galaxy scalelength Di we used tne radius containing 75% of 



the total (baryonic dark matter) mass of the system; 
- the parameter Tpjtenc that gives an estimate of the eff'ect of 
tidal forces exterted on the pair, per unit of time. 

We want to recall that, while for flybys this analysis relates 
the SFR at pericenter passage with quantities evaluated at the 
pericenter too, this is not the case of mergers, looking in this 
case for possible existing correlations between the SFR in the 
coalescence phase and physical quantities evaluated at the first 
pericenter passage. 

4.7.1 . Maximum SFR versus Rr, 



Fig. 21 shows the maximum SFR , as a function of the distance 
Rp at first pericenter passage, for flybys (lower panel) and merg- 
ers (upper one). 

It is not surprising to see that for flybys there is a tendency to 
have increasing star formation levels as the pericenter separation 
of galaxies in the pair diminishes. More interesting is the case of 
mergers, where a positive correlation (+0.7) between these two 
quantities is found, so that galaxies that at first passage are more 
distant are also those that suffer the most intense burst of star 
formation in the merging phase. 



- the tidal parameter Tp of the galaxy pair, defined a^ 



Tp = + Tp2, 



where 



Mr \R 



1,2 



(9) 



(10) 



quantifies the eff'ect of tidal forces at pericenter passage Rp, 
suff'ered by a galaxy of mass M/ and scalelength Di, due 



^ Note that in the following definition of Tp no dependence on the 
galaxy mass ratio is proposed. The sample being composed only by 
equal mass galaxies, this has no influence on the present study. A pos- 
sible formulation for unequal mass encounters could consist into intro- 
ducing a weight g = Mi /(Mi + Mcomp) in the expression for Tp^ i , i.e. 
Tp, i = logiQ[g{M,omplMi) X (Di/Rpfl 



4.7.2. Maximum SFR versus Vp 

No correlation is found between SFR maximum and the relative 
velocity of the two galaxies at the first pericenter passage for 
the merger sample, while a weak correlation (-1-0.34) exists for 
flybys (see Fig [22]). 



4.7.3. Maximum SFR versus l, 



As it is shown in Fig 23 the duration of the encounter, so, ulti- 
mately, the duration of the perturbation induced by the compan- 
ion galaxy, is a parameter that influences the peak of star for- 
mation (the absolute value of the correlation coefficient is about 
0.6 for both flybys and mergers). And, as previously found for 
Rp, the two samples of flybys and merger galaxies show a com- 
pletely difl'erent behaviour: indeed, while the maximum SFR in- 
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Fig. 21. Maximum star formation rate, relative to the isolated 
case, as a function of the first pericenter distance between the 
two galaxies , for flybys (lower panel) and mergers (upper 
panel). Difl'erent symbols correspond to difl'erent morphologies 
of interacting galaxies, as explained in the legends. The dashed 
lines represent the best linear least- square fits. Error bars repre- 
sent the standard error of the mean. 



Fig. 22. Maximum star formation rate, relative to the isolated 
case, as a function of the relative velocity at pericenter, for fly- 
bys (lower panel) and mergers (upper panel). Difl'erent symbols 
correspond to difl'erent morphologies of interacting galaxies, as 
explained in the legends. The dashed lines represent the best lin- 
ear least- square fits. Error bars represent the standard error of the 
mean. 



4.7.5. Maximum SFR versus Tpjte, 



creases as tenc diminishes for flybys, the opposite is the case for 
mergers, having a tendence to have higher bursts of star forma- 
tion at increasing tenc- This is exactly the same trend shown in 
Fig. 21 because of the fact that flybys with small pericentric 
distances must also have the highest relative velocities, and so, 
ultimately, small tenc- 



4.7 A. Maximum SFR versus Tp 

For mergers, a strong anti-correlation (-0.72) is found when plot- 
ting the maximum SFR versus the tidal parameter Tp, defined as 
the sum of the tidal parameters Tp exerted by the two galaxies. 
The clear trend in Fig. 24 is that galaxy pairs that suffer intense 
tidal effects at first pericenter passage tend to have lower star 
formation rates in the merging phase. The same behaviour is 
found when comparing the SEE at the merger epoch with Tp 
(see Fig. [25]). The opposite trend stands for flybys, i.e. at peri- 
center, galaxy pairs that suffers intense tidal effects tend to have 
also higher star formation rates. 



The same trend is found when comparing the SFR at the merger 
with T pi tenc ^ which gives an estimate of the perturbation due 
to tidal forces per unit time. With a correlation coefficient r = 
-0.55, it shows that the maximum SFR is somewhat affected by 
the tidal forces exerted by the companion galaxy per unit time. 
We want to recall that our star formation law is only density 
dependent, so that an increase in the SFR corresponds to an in- 
crease in the local gas density. In this sense, the anticorrelation 
found between the amplitude of the burst and the tidal forces 
exerted per unit of time can be understood if one considers that 
the strongest the forces, the greatest the quantity of galaxy mass 
(and so, ultimately, also of gas material) removed from the sys- 
tem, as effect of energy injection (see next Section for a deeper 
discussion). 

4.8. Why does the SFR in the merging phase depend on Tp ? 

The goal of this Section is to further investigate the relation 
found in Fig. 24 between the value of the SFR peak in the 
merging phase and effects of tidal forces on the pair, evaluated 
at the first pericenter passage. In all the previous sections, we 
have always refered our analysis to the couple of interacting 
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Fig. 23. Maximum star formation rate, relative to the isolated 
case, as a function of the characteristic time of the encounter 
tenc^ for flybys (lower panel) and mergers (upper panel). Diff'erent 
symbols correspond to diff'erent morphologies of interacting 
galaxies, as explained in the legends. The dashed lines represent 
the best linear least-square fits. Error bars represent the standard 
error of the mean. 



Fig. 24. Maximum star formation rate, relative to the isolated 
case, as a function of the tidal parameter Tp (see text), for fly- 
bys (lower panel) and mergers (upper panel). Difl'erent symbols 
correspond to difl'erent morphologies of interacting galaxies, as 
explained in the legends. The dashed lines represent the best lin- 
ear least- square fits. Error bars represent the standard error of the 
mean. 



galaxies: for example, the SFR is the total SFR, to which both 
galaxies of the pair contribute; the tidal parameter Tp has been 
evaluated following Eq.[TO| so summing the tidal eff'ects suff'ered 
by the two galaxies. In this Section, we will "decompose" our 
analysis, evaluating some physical quantities directly on each 
galaxy of the pair. This in order to have a more detailed picture 
of the role played by each component in the evolution of the 
interacting system. 

Being interested in explaining the relation found for mergers, 
here we restrict our analysis to the sample of merger galaxies, 
i.e. to all the galaxies involved in interactions that lead to 
mergers in 3 Gyr of evolution. 

First of all, for each galaxy in this sample, we evaluate the 
radius rjs containing 75% of the total mass (gas stars -hdark 
matter) of the system, just after the pericenter passage with the 
companion, this in order to quantify the expansion experienced 
by the galaxy, as an eff'ect of the encounter. In Fig 27 rjs is 
plotted versus the tidal parameter Tp^gai, which now is evaluated 
on the single galaxy and not on the pair. For simplicity, we 
have grouped galaxies into two subsamples: early-type systems, 
including gEO and gSa, and late-type ones, including gSb and 
gSd galaxies. 



Not surprisingly, the figure shows a strong correlation 
between these two quantities, indicating that the stronger the 
interaction at the pericenter passage, the greater the subsequent 
expansion of the outer parts of the system. The radius r^^ is 
computed including all the galaxy components, but it is likely 
that a similar relation holds for the gaseous component only, in 
the sense that the gas mass ejected into tidal tails grows as tidal 
eff'ects on the galaxy become stronger. This is shown in Fig [28^ 



where M, 



gas, out: 



the gas mass outside 20 kpc from the galaxy 



center, is plotted versus Tp^gai. M gas, out has been evaluated just 
after the pericenter passage and it has been normalized to the gas 
mass present, at the same time, in the disk of the corresponding 
galaxy evolving alone. 

To summarize, after the pericenter passage, the amount of 
gas in the galactic disk is depleted by two complementary ef- 
fects: 

1 . an increase in the star formation rate in the disk, which ac- 
celerates the gas consumption; 

2. the ejection of gas material from the disk into tidal tails. 

The role these two phenomena play in depleting the gas mass 
in the disk is shown in Figs. 29 and[30| 

The first figure shows Mgasjn, the gas mass inside 20 kpc from 
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Fig. 25. Maximum star formation efficiency, relative to the iso- 
lated case, as a function of the tidal parameter Tp (see text), for 
flybys (lower panel) and mergers (upper panel). Difl'erent sym- 
bols correspond to difl'erent morphologies of interacting galax- 
ies, as explained in the legends. The dashed lines represent the 
best linear least-square fits. Error bars represent the standard er- 
ror of the mean. 



the galaxy center, versus Tp^ 



p,gal 



. As in Fig. |28| Mgasjn has been 



evaluated just after the pericenter passage and it has been nor- 
malized to the gas mass present, at the same time, in the disk of 
the corresponding galaxy evolving alone. 
The second figure shows the dependence on Tp^gai of 



^gas- 



the integrated gas mass trasformed into star from the time ^ = 
to the time t = t(Rpj^ This quantity also has been normalized to 
the gas mass present, at the same time, in the disk of the corre- 
sponding galaxy evolving alone. 



Comparing these two Figures with Fig|28j it is evident that usu- 
ally the disk gas mass is decreased mainly because of the ejec- 
tion of a conspicuous gas fraction far out from the disk, while the 
encounter-enhanced star formation contributes only modestly to 
clear it out. Obviously, part of the gas mass lost in the tails, or 
acquired by the companion galaxy, can be reaccreted in the last 
stages of the merging event, thus contributing to the fuel avail- 
able in the burst phase. But this is not the case for all the gas mass 
lost from the disk. In this sense. Fig 28 contributes to explain the 
trend found in Fig|24](top panel). 



Fig. 26. Maximum star formation rate, relative to the isolated 
case, as a function of Tp/tenc (see text), for flybys (lower panel) 
and mergers (upper panel). Difl'erent symbols correspond to dif- 
ferent morphologies of interacting galaxies, as explained in the 
legends. The dashed lines represent the best linear least-square 
fits. Error bars represent the standard error of the mean. 



4.9. Toward a formulation for the SFE 

For practical use, it is convenient to have a fitting formula for 
the estimation of the SFE, even if limited to planar encounters. 
The aim of this Section is to propose a possible formulation, for 
mergers and flybys, as a function of: 1) the type of the encounter 
(E-E, E-L, L-L); 2) the orbital spin, if parallel or antiparallel to 
the galaxies angular momentum (i.e. direct or retrograde orbit); 
3) the pericenter distance R2b of the Keplerian orbit of two points 
with masses equal to the masses of the interacting systems. 

In other terms, we are looking for a function 
/(morph, spin, R2b) so that 



logioSFEmax = /(morph, spin, R2b). 



(11) 



As in the previous Sections, for mergers and flybys, the max- 
imum SFE is relative to that of the corresponding galaxies 
evolved in isolation. 



It is the time corresponding to the pericenter passage. 



In Section 4.7.4 we have shown the strong correlation exist- 
ing between the maximum SFE and the tidal parameter Tp at the 
pericenter, giving a global fitting formula, which takes into ac- 
count all the various morphological galaxy types. In general, this 
strong correlation between Tp and SFE still holds when group- 
ing galaxies in diff'erent classes depending on their morphologi- 
cal type and the orbital spin, direct or retrograde. This suggests 



16 



R Di Matteo et al.: Star formation efficiency in galaxy interactions and mergers: a statistical study 







1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
; ^ ^ early-type 


1 1 1 1 1 






- n n late-type 






o 

CO 


T 














o 

Oi 
















o 








T— 1 


, , 1 , , , , 1 , , , , 1 , , , , 


1 , , , , " 



-0,5 0,5 1 1,5 



Fig. 27. Radius containing the 75% of the total (baryons+dark 
matter) mass of the galaxy versus the galaxy tidal parameter 
Tp,gai- Both quantities have been evaluated just after the first 
pericenter passage. Diff'erent symbols correspond to early and 
late-type galaxies, as explained in the Figure. Red color refers 
to direct encounters, blue to retrograde ones. The analysis has 
been restricted to galaxies in the merger sample (see text). The 
dashed lines represent the best linear least- square fits. Error bars 
represent the standard error of the mean. 
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Fig. 29. Gas mass inside 20 kpc from the galaxy center, versus 
the galaxy tidal parameter Tp^gai. Mgasjn has been evaluated just 
after the pericenter passage and it has been normalized to the gas 
mass present, at the same time, in the disk of the correspond- 
ing galaxy evolving alone. Diff'erent symbols in the figure cor- 
respond to early and late-type galaxies, as explained. Red color 
refers to direct encounters, blue to retrograde ones. The dashed 
lines represent the best linear least- square fits. Error bars repre- 
sent the standard error of the mean. 
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Fig. 28. Gas mass outside 20 kpc from the galaxy center, ver- ^^8- ^^^^ * ^^^^^^ ^^e galaxy tidal parameter 



sus the galaxy tidal parameter Tp^gai. Mgas,out has been evaluated 
just after the pericenter passage and it has been normalized to the 
gas mass present, at the same time, in the disk of the correspond- 
ing galaxy evolving alone. Diff'erent symbols in the figure cor- 
respond to early and late-type galaxies, as explained. Red color 
refers to direct encounters, blue to retrograde ones. The analysis 
has been restricted to galaxies in the merger sample (see text). 
The dashed lines represent the best linear least- square fits. Error 
bars represent the standard error of the mean. 



gas- 



p,gal- 

represents the integrated gas mass transformed into star 



from the time ^ = to the time t = t{Rp). As in Figs 28 and 



29j this mass has been normalized to the gas mass present, at 
the same time, in the disk of the corresponding galaxy evolving 
alone. Diff'erent symbols in the figure correspond to early and 
late-type galaxies, as explained. Red color refers to direct en- 
counters, blue to retrograde ones. Error bars represent the stan- 
dard error of the mean. The analysis has been restricted to galax- 
ies in the merger sample (see text). 



to adopt a right-hand side in Eq[TT] containing a parameter that 
could approximate Tp. Note indeed that, for a practical use, we 
are looking for a formulation for SEE that could be expressed 
only in terms of the initial parameters of the simulated encounter 
and that does not depend on quantities that require an ad hoc 
simulation to be evaluated. 

We remind the reader that the tidal parameter Tp at pericenter 
has been defined in Eq|9]as the sum of the tidal parameters i 
and Tp2 relative to the two galaxies involved in the interaction, 
i.e. 



Rp being the pericenter distance and Di and D2 the galaxies 
scalelength. Rp was evaluated as the minimum distance between 
the two galaxy density centers and both Di and D2 were approx- 
imated with the radius containing 75% of the total galaxy mass 
and where evaluated at the pericenter passage. 
On the contrary, the suggested approximation for Tp will make 
use of all quantities yet available ab initio, i.e. 



Tp^ x = log 



10 



(14) 



Tp = logiQ 



log 



10 



M2 
Ml 



+ logiQ 



Ml 

M2 



D2 

Rn 



(12) 



(13) 



R2b being the pericenter distance of the unperturbed Keplerian 
orbit of two bodies, having masses equal to the two galaxy 
masses, and di and J2 being the radii containing 75% of the 
total mass of the unperturbed galaxies. 
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Table 4. Values for the A and B constants in Eq 15 for mergers 
and flybys. 



Mergers 



Flybys 





A 


B 


A 


B 


E-E ret 


-0.098 


1.557 


0.135 


0.653 


E-E dir 


-0.126 


0.826 


0.115 


0.724 


E-L ret 


-0.202 


1.089 


0.078 


0.482 


E-L dir 


-0.098 


0.779 


0.024 


0.543 


L-L ret 


-0.121 


0.954 


0.027 


0.461 


L-L dir 


-0.026 


0.621 


0.004 


0.467 



The suggested formulation for the SFE is the following 



logioSFE^ax = A log 



10 



R 



2b 



+ 5, 



(15) 



being A = A(morph, spin) and B = 5(morph, spin) constants 
depending on the morphology of the two galaxies involved in 
the interaction and on the orbital spin. Note that also the x term 
depends on galaxy morphologies, via di and d2. 

In Table |4j A and B are given, both for mergers and fly- 
bys, and for difl'erent classes (E-E, E-L, L-L, in retrograde or 
direct orbits). Their values have been obtained performing a lin- 
ear least-square fit between the x values and the corresponding 
values of logioSFEmax resulting from simulations. 
A comparison between the curves 



y = Alog 



10 



+ B 



(16) 



and the 'experimental' data for the SFE, for mergers and fly- 
bys, is shown in Figs [31] and [32] respectively. The proposed for- 
mulation for the SFE reproduces the data quite well, with the 
modulus of the relative errors quite uniformly distributed with a 
mean value that is always below 30%, with the only exception 
of Early-Late type direct mergers and Late-Late type retrograde 
flybys, where it is ^ 50%. 
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Fig. 31. SFE versus the parameter x defined in Eq 14 for merg- 
ers. Diff'erent panels refer to diff'erent classes for the encounters, 
depending on the morphology of the interacting galaxies and on 
the orbital spin. Error bars represent the standard error of the 
mean. The dashed line in each panel represents the proposed for- 
mulation for the SFE (see Eq[T5]). 
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4.10. Evolution in the (Lgas, ^sfr) plane 

In the previous sections, we have seen that the star formation 
rate in the merging phase anticorrelates with the intensity of 
the tidal interaction at first pericenter passage and we explained 
this in terms of gas removal from the galactic disk of the two 
galaxies, just after the encounter. In this picture, the next step is 
to verify if the SER during interactions and in the coalescence 
phase correlates with the quantity of gas material available in 
the galactic disk. 

Several observations of normal galaxies suggest that the 
SER on global scales is well represented by a Schmidt law of 
the form I^sfr = ^^gas^, ^gas being the gas density weighted 
inside the galactic optical radius and I^sfr the disk-averaged 
SER surface density. The best fitting slope N, obt ained with 
a conv entional least squares fit, is about 1.4 (see Kennicutt 
1998a). Even if with large variation for a given value of l^gas, 
these data show that a correlation between SER and gas density 
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for flybys. 



^ The sample being composed only by giant galaxies, in the following 
formula no dependence on the galaxy masses is proposed. We refer the 
reader to following papers for a more general formulation, that could 
take into account also the mass ratio of the interacting galaxies. 



Fig. 32. SFE versus the parameter x defined in Eq 
Diff'erent panels refer to diff'erent classes for the encounters, de- 
pending on the morphology of the interacting galaxies and on the 
orbital spin. Error bars represent the standard error of the mean. 
The dashed line in each panel represents the proposed formula- 
tion for the SFE, described (see Eq[T5]). 



exists, which is also largely independent of galaxy type. A 
similar correlation, with a comparable N slope, exists for 
starburst galaxies, if averaging the SER density and the gas 
density on the inner (1-2 kpc) galactic region, both the Ugas 
and the I^sfr values being in this case 1-4 orders of magnitude 
higher than those obtained for normal galaxies, thus suggesting 
that a large amount of gas is present in the central region to 
sustain the burst activity. 

Other authors found s ome variations in the value of the N slope: 
Wong & Blitz (2002]) studing a sample of seven molecule rich 
spiral galaxies found ~ 1.1-1.7, depending on the correction 



for extinction in Ha emission, in deriving the SER; Boissier et] 
[aL] p003 ) found - 2.0 for a sample of sixteen spiral galaxies; 
finally |Gao & Solomon] p004| ) suggested a star formation law 
with a power-law index of 1.0 in terms of dense molecular 
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Fig. 33. The global Schmidt law in galaxies. Upper left panel: (Zgas,^sFR) relation for interacting galaxies (grey empty squares). 
For comparison, also isolated galaxies are shown (stars). Upper right panel: (Zgas,^sFR) relation for merging galaxies. Different 
symbols in this and in the next two panels are explained in the legends. For comparison, the interacting galaxy sample is shown 
with grey empty squares. Lower left panel: circumuclear starburst. For comparison, interacting galaxies (grey empty squares) and 
mergers (black squares) are also shown. The dashed line represent a "global" Schmidt star formation law of the form I^sfr ^ 
l^gas^'^ ( |Kennicutt|l998a| , and it is shown for a comparison with the data. Lower righ panel: (Lgas,^sFR) relation for post-merger 
galaxies. For comparison, interacting galaxies (grey empty squares), mergers (black squares) and circumnuclear starburst (black 
empty circles) are also shown. 



gas content, studying a sample of 65 normal spirals, luminous 
infrared and ultraluminous infrared galaxies. 

In this framework, it is interesting to study if our sample 
obeys a similar global relation, provided that on local scales our 
star formation recipe is a Schmidt type. In particular, we are 
interested in understanding if interacting and starburst galaxies 
follow a "global" Schmidt law of the type I^sfr = ^^gas^ 



as found observationally for normal spiral galaxies ( [Kennicutt 
|1998a) and if, more generally, it is possible to trace a global 
evolution on the Ugas - fr plane, for pre-interaction, interact- 
ing, starburst galaxies and mergers. 

To compute Ugas and I^sfr, we evaluated the SFR and the 
gas amount inside a radius containing 85% of the total vis- 
ible (gas-hstars) galaxy mas^ and then we calculated the 
disk-density relative quantities. For interacting galaxies, we 



Note that for interacting and starburst galaxies, due to the fact that 
part of the mass populates tidal tails, thus spreading for several hundreds 
kpc outside the galactic center, we restricted our analysis to 85% of the 



evaluated the two quantities for each galaxy of the pair, while 
for starburst, circumnuclear starbursts and post-merger Ugas and 
the Us FR are evaluated on the resulting merger. 



The main findings of our analysis are summarized in Fig 33 
which shows diff'erent galaxy samples (interacting, starbursts, 
circumnuclear starbursts and mergers) in the plane (l^gas, fr)- 
In this plane the interacting galaxies (upper left panel) lie in a 
region which extends over 2 orders of magnitude both in llgas, 
from 0.1 to 30 M^pc-^, and in Y^sfr (lO""^ - OmU^yr-^kpc-^). 
No evident diff'erent behaviour is found for diff'erent morpho- 
logical type, i.e. for early and late-type systems. This sample 
includes galaxies at very diff'erent stages of interactions: there 
are systems that are well before the first close passage, hundreds 
of kpc distant from the companion (in this sense they can be 
considered as isolated) and systems close to the merging phase. 
We eliminated from this sample (as in the isolated sample 



total visible (gas+stars) mass located inside 20 kpc from the galaxy 
center. 
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shown in the same panel for comparison) the transient initial 
burst phase, corresponding to the emergence of density waves 
in the galactic disks (see Figs. [6] and [7]). Note that the best linear 
fit reported in Fig [33] for interacting galaxies would be steeper if 



the low ILgas - FR values in this sample were not considered. 

The starburst phase is represented in the upper right panel. 
The high correlation found for this sample (0.9) is remarkable. 
One can note that in this phase galaxies occupy an extended 
region in ^gas, and reach higher levels of SFR density, with 
respect to interacting pairs. Note also that mergers resulting 
from encounters between an elliptical and a spiral have, in 
general, higher gas surface densities and higher star formation 
rate density with respect to spiral-spiral mergers. The situation 
changes if one restricts the analysis of starburst galaxies to 
circumnuclear regions, in the inner kpc region of the galaxy. 
In this case (lower left panel) the nature of the starburst event 
clearly emerges: the gas circumnuclear density is, on average, at 
least 100 times higher than that of the overall disk, thus giving 
rise to a strongly enhanced SFR, whose circumnuclear density 
is far greater than those of interacting and starburst galaxies as 
a whole. In this case, late-type spirals encounters show higher 
SFR densities, if compared to mergers involving an early-type 
system. 

Once the starburst phase is over, the mergers evolve toward 
the lower left part of the (J^gas^^sFR) plane (as shown in the 
lower right panel in Fig 33 ), having at this point consumed most 
of their gas reservoirs. 

This overall trend in the location of normal and starburst galax- 
ies i n the (J^ggs^^sFR ) plane is confirmed also observationally 
(see Kennicutt]! 1998a , Fig. 9 left). Note that in this case highest 
values (for normal and merging galaxies) are reached both in 
the X and y axis, but it is just due to the fact that our simulations 
start with initial conditions typical of present day galaxies and 
then evolve for 3 Gyr in time, so that the amount of gas in the 
simulated systems is only initially comparable to that of galaxies 
in the local Universe. In particular, the position of isolated and 
interacting galaxies in Fig[33]top left is a reflect of gas amount. 



The present study on the evolution of isolated, interacting 
and merging galaxies in the (^gas.^sFR) plane extends and en- 
riches prior numerical works concerning star formation in iso- 
lated galaxies and in major interactions. Springel ( 2000), for ex- 
ample, studied the composite Kennicutt law for a small set of 
simulated isolated and merging galaxies, showing that his sam- 
ple remarkably reproduces the global Schmidt law over a large 
dynamic range. But the accordance of the data with the global 
Schmidt law is a natural consequence of the chosen star forma- 
tion and feedback recipes he adopted, being the free parameters 
chosen in order to satisfy Kennicutt 's findings. 
In an extensive study conducted to analyze how star formation 
depends on the adopted parametrization for feedback, |Cox et 
[aL] ( [2006| ) showed that isolated galaxies satisfy the empirical 
Kennicutt law, when including a gas density threshold for star 
formation. Indeed, when averaging star formation rate and gas 
surface density within azimuthal annulli, and adopting not too 
high feedback parameters, they found a good accordance with 
the results in Fig. 3 in ,Kennicutt| ( ^1998b| ). Then, studying the 
Kennicutt law for the merger sample, using an azimuthal aper- 
ture of radius 2 kpc, they confirmed that also the interacting sam- 
ple closely tracks the empirical star formation law. 
But both these works use a star formation and a feedback mod- 
elization designed to reproduce Kennicutt empirical law. In this 



sense, the present study investigates the ILgas - fr correlation, 
in a more general context. Indeed, starting from a star forma- 
tion recipe based on the local gas volume density (Eq[5]), it is not 
immediately clear that a global (on 10-20 kpc scales) relation, 
based on gas surface density, must hold. Evidently, we do not 
claim to reproduce the slope found for the empirical law, but, 
neverthless, it is striking to see how a global (i.e. on kpc scales) 
ILgas-^s FR relation is satisfied, for isolated, interacting, merging 
and post-merger galaxies. 



4.11. Gas Inflow 

It is well established that the interaction with a companion usu- 
ally destabilise the disk of a galaxy, which develops non axisym- 
metric structures (as spiral arms, bars), responsible of exerting 
torques on the gas material (Noguchi 1988, Barnes & Hernquist 
1996, Mihos & Hernquist 1996, Combes 2001). In this way, the 
gas loses angular momentum and falls i nto the inner galac tic re- 
gions, where a starburst takes place (see [Smith et al.|2007[ for an 
observational study of 35 interacting galaxy pairs, that show gas 
concentration in the inner regions, with a subsequent enhance- 
ment of the central star formation rate). 

In the previous section, we have shown indeed that the star for- 
mation enhancement in the merging phase is a strongly inhomo- 
geneous process, taking place mostly in the central kiloparsec. 
Here we want to describe in more details the gas dynamics dur- 
ing the diff'erent phases of the encounter. To do this, we studied 
the occurence and evolution of asymmetries in the disk stellar 
distribution and the response of the gas component to their emer- 
gence. To this aim, firstly, we Fourier- analyzed the disk surface 
density distribution, focusing in particular our attention to the 
amplitudes of the Fourier components m=l and m=2. Then we 
studied the evolution of the gas lagrangian radii containing 10%, 
25%, 50% and 75% of the total gas mass, in order to have a pic- 
ture of gas inflows into the central regions. Evidently, it is not 
possible to describe here the results of all the difl'erent encoun- 
ters simulated, so we will describe in detail only two of the cases 
analyzed. 

To Fourier-analyze the disk surface density distribution, we 
adopted the following procedure: 

- firstly, the angular momentum of the disk galaxy has been 
evaluated and the galactic plane has been rotated conse- 
quently, in order to have the galactic spin parallel to the z- 
axis of the reference frame; 

- then the x-y plane has been divided into an annular grid 
(exponentially spaced in radius and linearly in azimuth) be- 
tween the center C of the disk galaxy and \Qkpc. 

- the surface mass density distribution on the x-y plane has 
been evaluated, taking into account only the galactic old stel- 
lar component that lies between ±500 pc from the galactic 
disk. 

- once the density S(r, 6), being r the distance from the galaxy 
center and 6 the azimuth, is computed, for each annular ring, 
we have fitted E(r, 6) with the following function 



m 



1 -h \Aisin{6 + (pi) -h A2sin(26 -h ^2) + 



.... -\- A^sin(W -\- (p^)] 



(17) 



for each annu- 



obtaining a value of Aq, Ai,...,A8 and 0i, 
lar ring. 

finally, for each configuration, we have evaluated the aver- 
aged value < Al > and < A2 > of Ai and A2 respectively. 
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Fig. 34. Gas dynamics during the interaction and subsequent merger between an elliptical and a spiral Sbc galaxy. The first row refers 
to the two galaxies in direct orbit, the second row to galaxies in retrograde interaction. First column: relative distance between the 
two galaxy centers versus time. Second column: Evolution of the Ai strength, versus time (solid line). For comparison, the evolution 
of the Al asymmetry of the same galaxy evolving alone is shown (dotted line). Third column: Evolution of the A2 strength, versus 
time (solid line). For comparison, the evolution of the A2 asymmetry of the same galaxy evolving alone is shown (dotted line). 
Fourth column: Evolution with time of gas lagrangian radii containing 10% (black), 25% (red), 50% (blue) and 75% (green) of the 
gas mass are shown, in units of the corresponding gas radii of the isolated gSb galaxy. 



between 1 and 10 kpc and we adopted these values to quan- 
tify the "strength" of the m=l and m=2 asymmetries. 



The results of this analysis are shown in Fig 34 which refers 
to two encounters (a direct and a retrograde one) between an 
elliptical and a gSb galaxy. In the first column, the relative dis- 
tance between the two interacting system is shown, as a function 
of time, while the second and third column show, respectively, 
the evolution of m=l and m=2 asymmetries. It is evident from 
this plot that: 

- usually tidal encounters amplify both m=l and m=2 asym- 
metries; 

- at the first pericenter passage or during flybys, usually galax- 
ies in direct orbits develop A2 asymmetries more pronounced 
than those arising in retrograde encounters; 

- for mergers, a second amplification in the Ai and A2 asym- 
metries, relative to the isolated case, occurs when the two 
galaxies are in the final stage of coalescence. 



- in the phase between first pericenter passage and merger, the 
Al and A2 values can be either stronger or lower than those 
of the same galaxies evolved in isolation 

The emergence of strong asymmetries in the disk obviously 
influences the gas dynamics, as shown in the last column of 
Fig |34l w here the lagrangian radii for the gas component are 
showip] Indeed, after the first pericenter passage, as the disk 
galaxy developes transient m=l and m=2 asymmetries, a first 
gas inflow, involving about 25% of the gas mass takes place, ac- 
companied by an increase in the star formation efficiency. This 
inflow last about 200 Myr and it is followed by a re-expansion of 
the inner gas lagrangian radii, probably due to an enhancement 
in the gas kinetic energy from SNe explosions. But the most dra- 
matic and very rapid gas compression starts only when the two 
interacting galaxies are in the final phase of the merger: in this 
case the rapidly changing galactic potential causes the inflow of 



In fact here we refer to the gas+new stars component. 
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Fig. 35. Histograms of the maximum SFR (relative to the iso- 
lated case), for mergers. Flybys are shown for comparison in the 
small window inserted in the Figure. 



Fig. 37. Histogram of the relative ISFR, for mergers. Flybys 
are shown for comparison in the small window inserted in the 
Figure. 



dependence on several parameters. We have already pointed 
out (Section |4.5|) that interactions are not always a sufficient 
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Fig. 36. Histograms of duration of enhanced SFR, for the whole 
sample of interacting galaxies (mergers and flybys). Two thresh- 
olds are shown: relative SFR > 2 and relative SFR > 5 (the latter 
is shown in the small window inserted in the Figure). 



a high gas mass fraction (between 50% and 80% of the gas mass 
being involved) into the inner kiloparsec region, where a strong 
starburst takes place. 



5. Discussion: Are galaxy interactions always 
starburst triggers? 

In the previous Section, we have investigated in details the 
evolution of the star formation rate during flybys and mergers, 
showing the great variety of SFR evolutions occurring during 
galaxy interactions, and presenting a deep study of the SFR 



condition to convert high gas mass quantities into new stars. 
Here we want to deepen the discussion on this point, beacause 
of its potential impact also for observational studies. Evidently, 
we keep in mind that the statistical analysis presented in this 
paper contains inevitable limitations, that we will try to remove 
in following works. Even exploring a whole Hubble sequence 
(from early-tape to late-type giant galaxies), all the simulations 
have been performed using an unique dark matter model, for ex- 
ample. Having shown in Section [48] that tidal efl'ects are crucial 
in determining the SFR in the merging phase, it could be inter- 
esting to check, in subsequent works, how the overall analysis 
depends on the dark matter model adopted (varying the den- 
sity profiles, the limiting radii, the central concentrations, etc,...). 

The present work is also limited to coplanar encounters, with 
galaxies having null relative disk inclinations, and spins parallel 
to the orbital angular momentum. 

Anyway, the large sample (216) of interacting galaxies here pre- 
sented is yet sufficient to draw some preliminary conclusions 
about the interaction- starburst connection. 



Fig 35 shows the histogram of the maximum SFR for merg- 
ers (upper panel) and flybys (lower panel). As in the previous 
Section, this value is relative to the SFR of the galaxy evolved 
in isolation. The plot clearly puts in evidence that mergers are 
not always starburst triggers: indeed the number of mergers that 
produces star formation rates 10 times higher than those of iso- 
lated galaxy are about the 17% of the total merger sample, while 
half of the sample shows only a moderate enhancement in the 
SFR (less than 4 times the isolated case). 

High SFR are not only less frequent, but are also characterized 
by shortest duration times. This is shown in Fig 36 where the 
histogram of the burst duration is plotted, for two different SFR 
enhancement levels. About 36% of the whole interacting galax- 
ies samplq^ sustains SFRs two times higher than those of the 



i.e. including flybys and mergers together. 
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isolated counterparts for a time greater than 500 Myr (upper 
panel in Fig 36). But this duration significatively diminishes, 
for galaxies whose relative SFR is higher than 5. In this case, 
the maximum starburst duration is less than 600 Myr, and only 
13% of the whole galaxy sample sustains a SFR enchancement 
greater than 5, for more than 100 Myr (lower panel in Fig|36|). 
Finally galaxies that show SFRs 20 times higher than those of 
the isolated counterparts are able to substain it for 150 Myr, in 
the most favourable cases. 

Not only the starbust frequence is low, but also the integrated star 
formation rate, compared to the isolated case, seems to support 
the idea that interactions are not always sufficient to convert high 
gas mass quantities into new stars, as the histogram in Fig 37 
shows. Flybys and mergers can produce 2 times more stars than 
isolated galaxies, but about 50 % of all the mergers, and 76% of 
flybys, has an ISFR (see Section 



only 1.25 higher than that of iso 



6. Conclusions 



4.3.1 for its definition) which is 



ated galaxies. 



We have investigated the enhancement of the star formation rate 
in galaxy interactions, by numerical simulations, comparing 
the star formation properties of more than two hundred pairs 
of interacting and merging galaxies, with those of isolated 
galaxies. 

The present work partially confirms previous numerical inves- 
tigations (|Mihos & Hemquist||1994a|bl p^96| |Springel|[2QQQt 
Cox et al. 2006| ), in the sense that galaxy major interactions 
and mergers can trigger strong nuclear starburst, but it clearly 
puts in evidence that this is not always the case, i.e. mergers are 
not always starburst triggers and galaxy interactions are not 
a sufficient condition to convert high gas mass quantities into 
new stars. This mainly because strong tidal interactions at the 
first pericenter passage can remove a large amount of gas from 
the galaxy disks. This gas material, ejected into the tidal tails, is 
only partially re-acquired by the galaxies in the last phase of the 
merging event. 

We have shown that the star formation rate in the merging 
phase: 

1. does not depend mainly on the total amount of gas mass 
available just before the final coalescence phase, being this 
total amount of gas content the sum of that present in the tails 
and in the main body of the system. 

In turn, it is strongly anti-correlated with: 

2. the distance of the two galaxy centers at the first pericenter 
passage, in the sense that, on average, galaxies that suff'er too 
close passages produce also the lowest bursts of star forma- 
tion; 

3. the amplitude of the tidal forces at pericenter, i.e. pairs that 
suff'er less intense tidal actions at the first passage are able 
to preserve a great mass fraction in the disk, that consti- 
tutes the fuel for the nuclear starburst in the merging phase. 
Furthermore, 

4. the enhancement in the star formation in the merging phase 
depends on the galaxy spin, in the sense that, on average, 
galaxies in retrograde orbits are stronger starburst triggers 
than those involved in direct encounters. 

We have also analyzed the evolution in the (J^gas, fr) plane 
of interacting, mergers and post-merger galaxies, finding that, 
globally, the Kennicutt- Schmidt law is retrieved statistically for 
all the diff'erent stages of interaction. 



Finally, we have proposed a formulation for the SFE at the 
pericenter passage, for flybys, and in the coalescence phase for 
mergers. The general laws we derived depend only on the main 
parameters of the encounters, as the orbital spin, the pericenter 
separation and the galaxy dimensions. 

Obviously, it is still necessary to exploit a larger range of 
parameters, for example to get an insight into the dependence 
of the star formation efficiency on the orientation of the galaxy 
disks with respect to the orbital plane, or on the masses ratio of 
the galaxies in the pair. It would be interesting also to perform 
a subset of these simulations varying the star formation rule, i.e. 
from a density-dependent Schmidt type one, to a formulation 
that could take into account energy dissipation in shocks, as pro- 
posed by [Barnes ( 2004 >). In any case, in our opinion, this work 
can contribute, on one hand, to clarify the physical mechanisms 
behind the large interval of star formation enhancements found 
in observed interacting pairs, an to furnish, on the other hand, 
simple star formation formulations for theoretical modelling. 
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